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ABSTRACT 

Grains exposed to anisotropic radiation fields are subjected to forces due to the 
asymmetric photon-stimulated ejection of particles. These forces act in addition to the 
"radiation pressure" due to absorption and scattering. Here we model the forces due 
to photoelectron emission and the photodesorption of adatoms. The "photoelectric" 
force depends on the ambient conditions relevant to grain charging. We find that it is 
comparable to the radiation pressure when the grain potential is relatively low and the 
radiation spectrum is relatively hard. The calculation of the "photodesorption" force is 
highly uncertain, since the surface physics and chemsitry of grain materials are poorly 
understood at present. For our simple yet plausible model, the photodesorption force 
dominates the radiation pressure for grains larger than ~ 0.1/um exposed to starlight 
from OB stars. We find that the anisotropy of the interstellar radiation field is ~ 10% 
in the visible and ultraviolet. We estimate size-dependent drift speeds for grains in the 
cold and warm neutral media and find that micron-sized grains could potentially be 
moved across a diffuse cloud during its lifetime. 

Subject headings: dust 

1. Introduction 

Anisotropic radiation fields are common in the interstellar medium (ISM). Of course, the 
radiation is highly anisotropic at points located close to a star or cluster or at the edge of a 
molecular cloud. These include the astronomically important photodissociation regions (PDRs), 
the surface layers of molecular clouds, where incident UV dissociates molecules. The visible/UV 
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radiation field in the diffuse ISM is due primarily to starlight; since stars are unevenly distributed, 
it too is anisotropic. Dust grains in this radiation field are subjected to torques, if their shapes are 
somewhat asymmetric, and this may be resposible for the alignment of interstellar grains (Draine 
& Weingartner 1996, 1997), which is inferred from the observed polarization of starlight (Hiltner 
1949a, b; Hall 1949; Hall & Mikesell 1949). 

Grains exposed to anisotropic radiation are subjected to forces as well as to torques, again 
with potentially important consequences. If the "radiative" forces are large enough to result in 
substantial gas-grain drift, then this could lead to spatial variation in the dust-to-gas ratio. If 
the drift depends on grain size, then grain coagulation could result. Here we investigate two 
potentially significant forces, due to the recoil associated with electrons and adsorbed H atoms, 
which are ejected from the grain with some non-zero probability following photon absorption. 

Lafon (1990) modeled the net force due to anisotropic photoelectron emission for grains exposed 
to unidirectional radiation fields and found that, for some interstellar conditions, it substantially 
exceeds the radiation pressure, due to absorption and scattering. Lafon considered spherical grains 
and assumed that electrons are emitted only from the illuminated hemisphere. However, interstellar 
grain radii can be smaller than the photon attenuation length, so that electrons are also emitted 
from the non-illuminated hemisphere. In §3, we revisit the calculation of photoelectric recoil forces 
on dust grains, employing the prescription of Kerker & Wang (1982) for estimating the anisotropy 
in photoelectron momentum. 

In §4, we consider the force due to the anisotropic photodesorption of adsorbed H atoms. 
An adatom starts out in the gas phase, collides with and sticks to the grain surface, and can be 
removed when an absorbed photon breaks its bond to the surface. Of course, there are many 
other grain surface processes, besides adsorption and photodesorption, and these too could play 
important roles in determining the emission anisotropy. Since interstellar grains have not yet 
been well characterized, the details of their surface physics are unknown. We adopt a simple, yet 
plausible, model and evaluate the "photodesorption" force for this one case. 

Following the investigations of the photoelectric and photodesorption forces, we discuss, in §5, 
how surface processes can affect the drag force experienced by drifting grains. In §6, we estimate 
the anisotropy of the interstellar radiation field. In §7, we estimate drift speeds for grains in the 
diffuse ISM, tying together the results from the previous sections. In a separate paper (Weingartner 
& Draine 2001c), we will investigate gas-grain drift in PDRs. 



2. Radiation Pressure Force 

We generally normalize our results for the "photoelectric" and "photodesorption" forces to 
the radiation pressure force, in order to see whether or not they are significant. The force due to 
radiation pressure is 

F rad = 7ra 2 (Q pr )Au rad , (1) 
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where a is the grain radius (we assume spherical grains in all of our calculations) and cAu ra d is the 
net energy flux in the radiation field (Au ra( j = n ra( j, the total energy density, for a unidirectional 
radiation field and Au ra d = for an isotropic radiation field). The spectrum-averaged radiation 
pressure efficiency factor is given by 

, n , _ C" [Qabs + Qsca(l ~ (COS 9))] U v dv 

where Qabs^ra 2 is the absorption cross section, Q SC a,^o? is the scattering cross section, (cos 9) is the 
usual scattering asymmetry factor, and f max is the highest frequency in the radiation field. We will 
take /ifmax = 13.6 eV , because H I regions are highly opaque to H-ionizing radiation. We evaluate 
Qabsj Qsca> and (cos 9) using a Mie theory code derived from BHMIE (Bohren & Huffman 1983). 
We consider carbonaceous and silicate compositions and take optical properties from Draine & Lee 
(1984), Laor & Draine (1993), and Li & Draine (2001), as described by Weingartner & Draine 
(2001a). 

We consider two types of radiation spectrum. To characterize regions located near a single star, 
we adopt a blackbody spectrum (cut off at 13.6 eV ), with color temperature T c and dilution factor 
W, so that the energy density per unit frequency interval u v = AirW B U {T C ) / c. It is convenient 
to characterize the radiation intensity by G = ^rad/ u Hab' wnere u ^1d ls the energy density in the 
radiation field between 6eV and 13.6 eV and = 5.33 x 10~ 14 ergcm -3 is the Habing (1968) 
estimate of the starlight energy density in this range. 1 

For the diffuse ISM, we adopt the average interstellar radiation field (ISRF) spectrum in the 
solar neighborhood, as estimated by Mezger, Mathis, & Panagia (1982) and Mathis, Mezger, & 
Panagia (1983); see Weingartner & Draine (2001b) for a convenient representation of the ISRF. 
The total energy density in the ISRF is u = 8.64 x 10 -13 erg cm -3 , with u"^ d = 6.07 x 10 -14 erg cm~ 3 
in the 6-13.6 eV interval, or G = 1.13. 

In Figure 1 we display (Q pr ), as a function of grain size a, for the ISRF and blackbody spectra 
with various values of T c . We consider grains with a as large as l^xm, even though the classic MRN 
size distribution (Mathis, Rumpl, &; Nordsieck 1977) only extends up to a = 0.25/im. The size 
distributions of Weingartner &; Draine (2001a) contain non-negligible amounts of mass in grains 
with a ~ 1/xm, and we will find that the radiative forces are relatively large for such grains; thus, 
they could be important in some interstellar environments. We consider grains as small as a = 10A, 
equal to the value adopted by Weingartner & Draine (2001b) for the electron escape length, the 
typical distance that an excited electron travels within a grain before losing its energy. This choice 
is somewhat arbitrary, but this does not matter, since the radiative forces are inconsequential for 
grains with a ~ 10 A. 



1 For comparison, the interstellar radiation field estimated by Draine (1978) has u = 8.93 x 10 14 erg cm 3 between 
6 and 13.6 eV, or G = 1.68. 



-4 - 



3. "Photoelectric" Force 
3.1. Model 

We adopt the photoemission model of Weingartner & Draine (2001b). The photoemission rate 
and photoelectron energy both depend on the grain charge, so it is necessary to evaluate the charge 
distribution fz( a ) = the probability that the grain has charge Ze, where e is the proton charge. The 
charge distribution is set primarily by a balance between electron loss via photoemission and gain via 
accretion from the gas. Charging by this process is determined mainly by the parameter Gy/T/n e 
(T is the gas temperature and n e is the electron density), with an additional mild dependence on 
T (Bakes & Tielens 1994; Weingartner & Draine 2001b). 

The photoelectric force on a grain is given by 

F pc (a) = fz(a) [F pCiV (a, Z) + F pe , a (a, Z)] , (3) 
z 

where F pe ^(a, Z) is the force due to electrons emitted from the valence band. When Z < 0, the 
— Z "attached" electrons occupy energy levels above the valence band, if the latter is full in the 
neutral; F pe ^(a, Z) is the force due to the emission of these attached electrons. 

The contribution from the photoemission of valence electrons is 

Z)=vra 2 / dvAYQ ahs -^V2^~e dEf E (E)y/E S(Z,a,E) , (4) 

where the photoemission threshold energy hv pe t is the minimum energy that an absorbed photon 
must have for an electron to be ejected from the valence band, E m [ n (E max ) is the minimum 
(maximum) photoelectron kinetic energy (at infinity), and /e(E) is the photoelectron kinetic energy 
(at infinity) distribution; see Weingartner & Draine (2001b) for prescriptions for evaluating these 
quantities. The emission asymmetry factor A(hu, a) measures the asymmetry in the emission of 
photoelectrons over the grain surface, the recoil suppression factor S(Z, a, E) accounts for electron 
emission in directions other than the surface normal, and m e is the electron mass. 

When Z < 0, the contribution from photodetachment is given by 



F pe ,a,{a, Z) — 




dvAo vAt (y)—± L J'2m e ,{hv - hu pdt + E min )S(Z, a,E = hv - hv pAt + E min ) 



(5) 

where the photodetachment threshold energy /if p dt is the minimum energy that an absorbed photon 
must have for an attached electron to be emitted from the grain and cr p dt(^) is the photodetachment 
cross section. 

We follow the simple prescription of Kerker & Wang (1982) for estimating A. The probability 
of photoemission from any site on the surface is taken to be proportional to the electric intensity 
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|E| just below the surface at that point. Thus, 



-£sm9cose\B(e)\ 2 de 

Ainu, a) = ^— q 

J^sm^lE^)! 2 ^ 



(6) 



where 9 is the polar angle with respect to the direction of the incident radiation. The electric field 
E(#) is evaluated using Mie theory (Bohren & Huffman 1983). 2 In Figures 2 and 3 we display the 
asymmetry factor A(hu, a) for grains exposed to unidirectional radiation fields as a function of inci- 
dent photon energy for various grain sizes for carbonaceous and silicate composition, respectively. 
Generally, A(hv,a) equals the value for a unidirectional field times A« ra( j/u ra d. Curiously, A < 
for some values of (hv, a). 

We assume that the electrons emerge symmetrically with respect to the local surface normal, 
with a "cosinedaw" angular distribution (i.e. the emission rate at angle ip with respect to the surface 
normal oc sin ip cos ip dip) . For an uncharged grain, this would imply S = 2/3. When the grain is 
charged, the electron escapes on a hyperbolic trajectory, so that when it is at infinity its velocity 
vector makes an angle ipoo with respect to the surface normal which differs from the corresponding 
angle at the surface, ip. Taking this into account, we find 



S = 



J^ 2 dtp sin tp cos tp cos f/'oo 

J^ 2 dip sin ip cos ip 
2 —- b 6 2 (l + 6/2) 
3 VTT1) -2 + 2(1 + 6) 









In 













'1 + 6/2 - VTTb 
l + b/2 + y/TTB 



+ 



2^/TTb 



b = 



(Z + l)e 2 
aE 



(7) 



(8) 



where Ze is the grain charge prior to photoelectron emission. For — 1 < 6 < 10 4 , the simpler 
expression 

1.107 

~ 0.46 69 + (6 + 1.421) - 5043 ^ ' 

is accurate to within 0.6%; we adopt this expression in our computations. When an electron tunnels 
across the Coulomb barrier of a negatively charged grain, 6 < — 1 and the above classical derivation 
fails, since the grain surface is a classically forbidden location. We simply take S = 1 when 6 < — 1. 



3.2. Computational Results 

The photoelectric force depends on the grain charge distribution, and hence on the parameters 
G\[T jn e and T. In Figures 4 and 5, F pc is plotted for carbonaceous and silicate grains, respectively, 
for T = 10 2 and 10 3 K and three values of G\fTjn e , ranging from 10 3 to 10 5 K 1 / 2 cm . For these 



2 We adopt dielectric functions as in Weingartner & Draine (2001b). 
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cases, we adopt a blackbody spectrum with T c = 3 x 10 4 K, cut off at 13.6 eV . We also provide 
results for conditions appropriate for the cold neutral medium (ISRF, T = 100 K, n e = 0.03 cm -3 , 
GVT/n e = 380) and the warm neutral medium (ISRF, T = 6000 K, n e = 0.03 cm- 3 , GVT/n e = 
2900). 

Note that for a given value of GVT/n e , the results do not change much as the temperature 
increases from 10 2 to 10 3 K. Also, F pe /F Ta d decreases for higher Gy/T/n e because the higher 
ionization potentials quench the photoemission. 

For the truncated blackbody spectrum, there is no systematic increase or decrease of F pe /F ia d 
over the full range of grain sizes, because the increase in the anisotropy of emitted electrons with 
grain size roughly compensates for the decrease in photoelectric yield. The sharp minimum near 
300 A in Figure 4 results from a local plateau in the variation of the anisotropy with grain size, 
which lies near the maximum of -F r ad- For the ISRF, the photoelectric force becomes less important, 
compared with the radiation pressure, as the grain size increases beyond ~200A. As seen in Figure 
1, for the relatively soft ISRF, the radiation pressure increases rapidly as a increases from 10 A to 
1/zm. The total energy density between and 13.6 eV is about an order of magnitude greater for 
the ISRF than for a blackbody with T c = 3 x 10 4 K and G = 1. For a = 10 A, F rad is comparable 
for the two spectra, but is about an order of magnitude larger for the ISRF when a = 1/xm. Since 
G = 1.13 for the ISRF, F pc is comparable for the two spectra. 

In Figure 6 we show the effect of varying T c over a range from 2 x 10 4 to 5 x 10 4 K, characteristic 
of hot stars, using carbonaceous grains and GVT/n e = 10 3 K 1 / 2 cm , T = 100 K as an example. 

4. "Photodesorption" Force 

4.1. Introduction and Model 

In the process of adsorption, gas-phase atoms or molecules collide with and stick to a solid 
surface, forming either a weak van der Waals bond (physisorption) or a stronger chemical bond 
(chemisorption). Adsorbed atoms can be removed from interstellar grains by various desorption 
processes: 1) thermal desorption, in which a thermal fluctuation in the grain breaks the adsorp- 
tion bond, 2) photodesorption, in which an absorbed photon breaks the bond, and 3) exothermic 
chemical reaction between two adatoms. If photodesorption is important, we expect a non-zero net 
recoil force in an anisotropic radiation field. 

The details of adsorption on grain surfaces are poorly known at present. In particular, it 
is unclear whether adsorbed hydrogen atoms can diffuse across the surface, either by thermal 
barrier hopping or quantum mechanical tunneling. In their classic study of H2 formation on grains, 
Hollenbach & Salpeter (1971) assumed that H physisorbs with binding energy ~ 0.03 eV , implying 
high mobility. Physisorption occurs on ice surfaces, so that this description is likely valid deep 
inside molecular clouds, where grains are covered with ice mantles. However, grain temperatures in 
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PDRs are high enough that any ice mantles have sublimed, and ice features (e.g. at 3.1/itm for H2O 
ice) are not observed in the extinction curves for the diffuse ISM. Tielens & Allamandola (1987) 
discuss experimental evidence that H chemisorbs on graphite and amorphous silicate surfaces, with 
binding energy in the 1-2 eV range. Also, theoretical studies (Mendoza & Ruette 1989; Fromherz, 
Mendoza, & Ruette 1993) of H adsorption on graphite find chemisorption with binding energies in 
this range. For such large binding energies, thermal diffusion is negligible. The tunneling rate is 
approximately given by the oscillation frequency of the atom in the chemisorption potential times 
the barrier penetration factor, so that 

~2d 



Tt ~ 10 13 s exp 



h (2mE b f 2 



(10) 



where Tt is the tunneling timescale, d is the distance between binding sites, m is the mass of the 
adatom, and is the binding energy. The argument of the exponential is a large number; if (Ef,, 
d) = (1 eV , 1 A), then Tt ~ 13 days for H adatoms. If d is increased to 2 A, then T t ~ 4 x 10 17 years. 
We assume that no diffusion occurs, primarily because this case is easier to model. Once a binding 
site is occupied by a chemisorbed H atom, another H can physisorb on top of the first (Tielens & 
Allamandola 1987). An atom in such a physisorption state can easily migrate to neighboring sites; 
also, if it can overcome an activation energy, it can react with the underlying H atom. 

In our model, we assume that each surface site is a chemisorption site, from which adatoms 
cannot diffuse. For simplicity, we also ignore gas-phase species other than H. When an H atom 
from the gas arrives at an empty site, it sticks with probability sn; otherwise it is assumed to 
undergo perfectly elastic specular reflection. Since the binding energy in the chemisorption site 
greatly exceeds the kinetic energies of the atoms in the gas, we take sh = 1- Also, we make 
the simplifying assumption that an H atom arriving at an occupied site either combines with the 
resident chemisorbed atom, with probability s c , or else undergoes elastic specular reflection. This 
approximation is most reasonable for high gas temperatures, when the kinetic energies of the gas- 
phase atoms exceed the physisorption energy and reaction activation barrier. We also assume that 
H2 molecules leave the grain immediately upon formation. 

To calculate the net force on a grain we must evaluate, as a function of location on the 
grain surface, the average rate at which momentum is transferred to the grain per surface site. 
Momentum transfer occurs when an atom is photodesorbed, a newly-formed molecule leaves the 
grain, a gas-phase atom reflects off of the grain, or a gas-phase atom sticks to the grain (either at 
an unoccupied site, where it remains until it is removed by photodesorption or molecule formation, 
or at an occupied site, where it immediately forms a molecule, which leaves the grain). 

We denote the rate at which gas phase H atoms arrive at a site by i? arr . The rate at which 
incident H sticks per empty site is given by SH-^arr and the rate at which H2 is formed per occupied 
site is given by s c -R arr . When a grain drifts with respect to the gas the arrival rate i? arr varies 
over the grain surface; of course, this results in the drag force. Baines, Williams, k, Asebiomo 
(1965) calculated the drag, assuming that all arriving atoms undergo specular reflection. The 
surface processes considered here substantially complicate the drag calcuation and the variable 
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i? arr complicates the net force due to photodesorptions. We will call the net force on a non-drifting 
grain, due to all of the processes mentioned above, the "photodesorption" force. In §5 we will 
discuss how the drag force is modified when the surface processes are taken into account. 



4.2. Photodesorption Rate 

The rate at which H is photodesorbed per occupied site is denoted by R p d = Rp d ^(9'), where 
i?p d is the photodesorption rate averaged over the grain surface (« the incident photon flux density 
times the photodesorption cross section) and ty(0') describes the variation in photodesorption rate 
as a function of the angle 9' with respect to the direction S of net energy flux. 

The photodesorption cross sections are not known. There is a substantial literature on ex- 
perimentally determined cross sections for simple molecules on a variety of substrates, which is 
summarized in the review by Franchy (1998). For example, for NO and photon energies ~ 6eV or 
less, the cross sections range from ~ 10 -22 cm 2 to 2 x 10~ 17 cm 2 . The cross section increases with 
progression through the following sequence of substrates: platinum metals, noble metals, semicon- 
ductors, oxides. The differences are due to variation in the efficiency with which the substrates 
quench the excited states. 

Hellsing et al. (1997) studied the photodesorption of K from graphite for photon energies 
< 5 eV and found a maximum cross section of « 2 x 1CT 20 cm 2 . However, for these photon energies 
the desorption is substrate-induced. That is, a bulk excitation (e.g. a photoelectron) travels to the 
adsorbed atom and excites it. 

Photons with adequately high energy can directly excite the adsorbed atoms. We estimate 
the photodesorption cross sections due to this direct process, for H on graphite and silicates, by 
considering the photodissociation of the CH and OH molecules. Photodissociation cross sections 
for CH and OH have been calculated by van Dishoeck (1987) and van Dishoeck & Dalgarno (1984), 
respectively. The cross sections reach values as high as a few xlCP 17 cm 2 , and are appreciable for 
photon energies > 6eV. Roberge et al. (1991) found photodissociation rates of 8.6 x 10 _10 s _1 
and 3.5 x 10 _10 s _1 for CH and OH, respectively, for the Draine (1978) radiation field. In the 
absence of more relevant experimental evidence, we simply set the photodesorption rates equal to 
the photodissociation rates. Also, in our crude approximation, we ignore variations in the shape of 
the radiation spectrum, and simply take the photodesorption rate proportional to u^ d , the total 
energy density between 6 and 13.6 eV. Since G = 1.68 for the Draine radiation field, we take 
i?p d = 5 x 10 _10 Gs _1 for carbonaceous grains and R° pd = 2 x ltr 10 ^- 1 for silicates. 

To evaluate *$>(9') we assume that, for any point on the grain surface, R p d is proportional to 
the electric intensity |E e | 2 just above the surface at that point. Thus, 

2 !^f h /h du Uv \E e {9')\ 2 
Jo d ^ sm ^J6 C v/h dvu u \E e (£)\ 



-9- 



In Figures 7 and 8, we display ^(6 r ) for carbonaceous and silicate grains, respectively, for a unidi- 
rectional radiation field and several values of grain radius a. Here we assume a blackbody spectrum 
(cut off at 13.6 eV) with T c = 3 x 10 4 K. For a < 100 A, #(0') looks like the 100 A curve, with 
greater symmetry for smaller grains. Generally, 

^ = Kad ~ A^ rad ) + AlZ rad tt°(fl') 
^rad 

where ^°(6 r ) is the value for a unidirectional field. 



4.3. Surface Site Occupation Fraction 

The momentum transfer rates depend on the occupation fraction /, i.e. the fraction of the time 
that a site is occupied by an H atom. If the grain's surface is fixed with respect to the radiation, 
then / is a function of the angle 6' with respect to S. However, the polarization of starlight implies 
that grains spin, with their spin axes statistically aligned (presumably with the magnetic field). 
Purcell (1979) noted that systematic torques likely act on grains, resulting in rotational speeds 
orders of magnitude greater than what would be expected from equipartition arguments. The spin 
period is much shorter than the timescales on which the surface processes occur, so that / is a 
function of the angle with respect to the spin axis u. 

We evaluate /(0) by balancing the the rate at which H is removed from an occupied site by 
photodesorption and H2 formation against the rate at which H is added to an empty site when 
arriving gas-phase atoms stick: 

/ (R pd + S C R & „) = (1 - /) SH^arr , (13) 

where bars denote averaging over the spin. This yields 

f(e) = s n [su + s c + ay(9)]~ 1 . (14) 

The parameter a is defined by 

"-f 1 > (15) 
where the H arrival rate for a non-drifting grain is given by 

RL = ^ H / 2 (7r/r 1/2 ; (is) 



«h is the H number density, I is the area of a surface site, and (3 = m p /2kT. Of course, 

1 f 2n 

where 4> is the azimuthal angle about Cj. 



(17) 
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Lazarian & Draine (1999a,b) have recently suggested that grains sometimes undergo rapid 
"flipping" of the principal axis of largest moment of inertia (fixed in body coordinates), alternating 
between being parallel or antiparallel to the angular momentum vector. We consider the limit 
where flips occur instantaneously, so that 

f(0) = sn{su + s c + ^[^(e) + ^(7r-e)}y 1 . (18) 

Although rapid flipping does not seem likely for suprathermally rotating grains with a > O.l^m 
(because of the large energy barrier that must be overcome in order to flip), the calculation of 
the photodesorption force in this case is instructive. Thus, we will consider both flipping and 
non-flipping grains. 



4.4. Evaluation of the Photodesorption Force 

We assume that the direction of net energy flux S lies at angle # ra d with respect to the grain's 
spin axis a> 3 . We adopt a coordinate system with the z-axis lying along u and the x-axis oriented 
so that S lies in the x-z plane. The z-axis serves as the polar axis for spherical coordinates, with the 
x-axis the reference for azimuth; thus the spherical coordinates for S are (# ra d,0). The "radiation" 
angle 9' is given by 

cos 9' = sin 6* ra d sin 9 cos 4> + cos ra d cos 9 . (19) 

We assume that both the photodesorbed H atoms and newly-formed H2 molecules emerge 
symmetrically with respect to the local surface normal, with a "cosine-law" angular distribution, 
as adopted for photoelectrons. The momenta of the outgoing atoms and molecules are denoted by pu 
and pn 2 respectively. We ignore polarization of the outgoing atoms and molecules when the grain 
has non-zero charge. The momentum transfer, per surface site, due to each of the contributing 
processes is directed radially inwards, with rates as summarized in Table 1. The factors of 2/3 
appearing in Table 1 are due to the cosine-law angular distribution. For reflections and sticking, 
the term in brackets is the probability that the process will occur; we include the absorption of 
the gas-phase momentum of an arriving atom which will immediately undergo H2 formation in 
the "sticking" category. The probability that the component of a gas H atom's velocity along the 
surface normal lies between v and v + dv is given by P(v)dv = (/3/vr) 1 / 2 exp(- (5v 2 )dv . 



Integrating the momentum transfer rates over all surface sites yields the photodesorption force: 
" ! - -z / d9 sin^ cos 9f(9) [C + ¥(0)] -x / d9sm 2 9f(9)— dxj> cos <t>V{0') , (20) 

Jo JO 27T Jo 



77ra 2 



where 



3 We consider # rad to be constant; i.e. we ignore the likely precession of u> about the magnetic field. 
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and 



C=- s c 



a 




(22) 



Since the H arrival rate does not depend on surface location for a non-drifting grain, the 
contributions due to H2 formation, reflections, and sticking depend only on f{9) and are thus 
directed along z = u; they are all accounted for in the constant C. For rapidly flipping grains, 
f(6) = f(TT — 6), so that only photodesorption contributes to the momentum transfer; this is also 
true for non-flipping grains when 6* r ad = 90°. In these cases, for given {7, sr, s c , radiation field}, 
the maximum force occurs when «h — > 00, so that a — ► 0, f(0) — > constant, and 



When /(#) 7^ f(ir — 9), the magnitude of F pc j can be much greater than F^, depending on the 
value of C. For example, when a <C 1 and aC S> 1, -F p d ~ aCi^/(sH + s c ), for non-flipping 
grains. However, for realistic parameter choices and interstellar conditions, it is usually the case 
that F p( j < F^. Also note that, for non-flipping grains with 6> ra d = 0, F p( j = when aC = sh + s c . 
In this case, the contributions from the various momentum transfer processes exactly balance, 
regardless of the functional form of ^f(9'); of course, such a parameter conspiracy is unlikely for 
real grains in the ISM. 



In Figure 9 (lower panel), we display the photodesorption asymmetry integral that appears in 
equation (23) for F^, for carbonaceous and silicate grains. This integral takes nearly the same 
values for the ISRF as for a blackbody (cut off at 13.6 eV) with T c = 3 x 10 4 K. We found in 
Figures 2 and 3 that, for some grain sizes and photon energies, the electric intensity just below the 
grain surface is unexpectedly concentrated in the hemisphere facing away from the radiation. The 
same effect occurs with the external electric intensity, and leads to a "backwards" photodesorption 
force, for a range of grain sizes. 

In the top panel of Figure 9, we plot F^/ F T&( \ for our canonical parameter values, given in 

Table 2. The surface site area of / 2 = 6 A 2 would result if a site is a hexagon in a graphite sheet with 
C-C spacings of rs 1.5 A. We assume that photodesorbed H atoms carry kinetic energy En = 2 eV 
and newly-formed molecules carry kinetic energy Eu 2 = 1 eV , implying pu = pn 2 (note that pn 2 
is not relevant for the calculation of though). For large grains the photodesorption force 
dominates the radiation pressure for the blackbody spectrum, but is not nearly so important for 



In the upper panel of Figure 10 we display F p d/F^ for non-flipping carbonaceous grains with 




(23) 



4.5. Computational Results 



the ISRF. 
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a = 0.3/um, for a few values of a. We also assume that sh = s c = 1, pn 2 = Pk, and a blackbody 
radiation spectrum (cut off at 13.6 eV) with T c = 3 x 10 4 K. In the lower panel of Figure 10 we 
display cos#f, where 9{ is the angle between F p( j and S. The results for other grain sizes for which 
F p< i is significant, for silicates, and for illumination by the ISRF are very similar. 

As # ra d increases, H atoms which stick to sites on the "dark" side of the grain have greater 
opportunity to photodesorb from the side facing the radiation, as they are transported there by 
the rotation; thus the photodesorption force increases with $ rac [. For flipping grains, the occupation 
fraction / depends only mildly on 6 (Figure 11), so that F p( j ~ and cos#f ps 1 for all values of 

#rad- 

In Figure 12 we plot F p ^/F^ and cos#f for s c = 0.1 and various combinations of a and 
(kT/En) 1 / 2 ; otherwise the parameters are the same as for Figure 10. When a = (kT/E n ) 1 / 2 = 0.01, 
cos Of is very nearly 1 for all values of # ra d, much more so than when s c = 1 and a = 0.01. Although 
f(9) is somewhat more constant for the case with sji = 1, C is also substantially larger, and it 
accounts for the bulk of the component perpendicular to S. When (kT/En) 1 / 2 = 1, C is large 
enough that cos#f « — 1 for 6> ra d = 0. 



5. Drag Forces 

Baines et al. (1965) calculated the drag force on a spherical particle drifting through a di- 
lute gas, assuming that colliding atoms reflect elastically. Draine & Salpeter (1979) obtained an 
analytical expression which accurately reproduces the Baines et al. numerical results, and this is 
widely employed in studies of drifting interstellar grains. For subsonic drift speeds, the drag due 
to collisions with neutral atoms of element j, in the absence of surface processes, is given by 



FL S , 3 ^™ 2 °-n jVd ^ Ay , (24) 



' 2m p kT \ l l 2 1/2 

'drag,j - 3<^d ^ - j 

where rij and Aj are the number density and mass number of element j, respectively. 

The Coulomb drag, due to interactions with distant ions, is often important in interstellar 
applications. For subsonic drift speeds, the Coulomb drag is given by 



1 f m \ ' ( 

^drag.Coul ~ Fdrag,H ) E ^ [ 



rrij \ „ / Ze 2 \ 2 , 

In 



akT 



3(kTf/ 2 



2e 3 \Z\(TTxn K 



1/2 



(25) 



(Draine & Salpeter 1979). The sum is over ions i, and x = ^ x«. Since the grain charge fluctuates 
on a time scale short compared to the drag time, we average over charge states. The dominant 
ions in neutral gas are C + and H + . We take the gas-phase C abundance to be C/H = 1.4 x 10~ 4 



4 For our assumed carbonaceous grain parameter values, G/nn w 4.4 x 10 3 aVTK '' 2 cm 3 , 
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(Cardelli et al. 1996; Sofia et al. 1997). Thus, when calculating Fdrag.Couh we assume that xq = x, 
x H ~ when x < 1.4 x 10" 4 , and x c = 1.4 x 10" 4 , x H = x - 1.4 x 10" 4 when x > 1.4 x 10" 4 . 

Next, we consider how the surface processes considered in the previous section, namely H 2 
formation and photodesorption, affect the drag force Fd ra g,H due to collisions with atomic H. 

For a drifting grain, the arrival rate i? arr depends on the angle 9" with respect to the grain's 
velocity u: 

R arr (0") = Ri r D(e") , (26) 



where 



D{6") = exp i-y 2 ) + V^V 1 + 4= f V dt exp (-t 2 ) 

I Jo 

with y = cos 9" ' . The expansion 



(27) 



D{e") « l + v^y + y 2 - 



(28) 



is accurate to within 1% for —0.653 < y < 1.05. 



The momentum transfer rates due to reflections and sticking are modified from the non-drifting 
case by P(v)v 2 — ► P(v)(v + ucosO") 2 in Table 1 (and the lower integrand is now —ucos9"). The 
rates obtained for non-drifting grains must be multiplied by the factor 2E(9")/^, where 



E{9") = yexp(-y 2 ) + (l + 2y 2 



exp(— t 2 )dt 



^- + 2y + vW + ^y 3 - ^y 5 



(29) 



The expansion in equation (29) is accurate to within 1% for —0.749 < y < 1.43. 



Photodesorption substantially complicates the drag, since the occupation fraction is no longer 
constant, but is given by 



f(0) = SK 



su + s c + a 



D{9) 



(30) 



for non-flipping grains. Of course, for flipping grains ^>(9)/D(9) is replaced by [*(#) + ^(vr — 
0)\/[D{9) + D(tt — 9)] in equation (30). The drag force can be found by evaluating the total force 
resulting from collisions of H with the grain (i.e. F p d + Fdrag,H~) and subtracting F p d; the total 
force is given by 



a 2 r f 2w 

Fpd + F drag ,H = J q d9 sin 9 # 

^R° pd p H f(eM0') + \Rl„ Scm J{9)D{9") + [2 - s u + f(9)(s H - s c )} ^R° aTT E(0" 

[z cos 9 + x sin 9 cos <p + y sin 9 sin 4>] 



(31) 
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To first order in u, equation (31) yields, for nipping grains, 



* drag,H 
^drag.H 



2 - s H . 3 
u 



SH - S c + 



ScPH 2 



7T 



\ 1/2' 



2m p kT / 



d6 sin 0/ o (0) [sin 2 0u - (1 - 3 cos 2 0) cos u u>] , 

(32) 



where fo(0) is the occupation fraction that would hold for the same grain if u = 0. For non-flipping 
grains, the following term must be added to the right-hand side of equation (32): 



where 



and 



a cos U 
Ji(0) = z cos 



d0sin0cos0/ 2 (0)^(0) [Ji(0) + J 2 (0)] 



3-7T 



7T 



2m p kT 



1/2 



J 2 (0) = q^ h 



2m p kT 



1/2 



z cos 0*(0) + x sin 0— / # cos <£*(0') 
27r y 



(33) 



(34) 



(35) 



The above approximation works best for small a and a. For the range of conditions T 6 (50, 6000) K 
and y^u < 0.1, it is valid to within 1% when a < O.l^m, a < 1; 3% when a < 1/jm, a < 1; 10% 
when a < 0.1/im, a < 10; and 30% when a < lfim, a < 10. The direction of Fd ra g,H can differ 
substantially from — u when a > 1, depending on the relative orientations of u, a>, and S. 



6. Starlight Anisotropy 

Here we estimate the radiation field anisotropy in the diffuse ISM, by considering the solar 
neighborhood to be typical and adding up the flux density over the entire sky. Most of the flux in the 
visible and UV comes directly from stars; starlight scattered off of dust grains (the diffuse galactic 
light, DGL) makes a non-negligible additional contribution (see, e.g., Witt 1989). Ideally, we would 
like to have a catalog of stars with accurate photometry and complete down to the magnitude 
beyond which fainter stars do not significantly contribute, along with accurate photometry of the 
DGL. 

Although this ideal is not currently available, there are a few catalogs from which useful 
estimates can be made. The Skymap Star Catalog, Version 3.7 (Slater & Hashmall 1992) is claimed 
to be complete down to magnitude 9 in the B and V bands, and the Tycho Catalog (ESA 1997) is 
claimed to be 99.9% complete down to V ~ 10.5. However, the photometry is missing for thousands 
of stars in the Tycho Catalog, including many bright stars. In Figure 13 we plot the integrated 
energy densities per unit frequency u u , normalized to the ISRF, as a function of the limiting stellar 
magnitude. The Skymap results are for the Johnson B and V bands (Bj and Vj respectively); the 
Tycho bands Bj> and Vr closely resemble the corresponding Johnson bands. The Tycho catalog 
also gives indirectly-derived values for Vj for all stars, including those for which Tycho photometry 
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is missing. The Bt and Vr energy densities are systematically lower than the Bj and Vj energy 
densities due to omission of stars lacking Tycho photometry. The Vj fluxes from the two catalogs 
agree very well, except at the limit of the Skymap catalog. Note that at magnitude 10 the curves 
are still concave up, with the integrated starlight accounting for less than half of the ISRF. 



The normalized starlight dipole moment is given by 



u^ISRF) 



(36) 



where u u (s) is the energy density per frequency interval due to star s and n(s) is the direction to 
star s. In Figure 14 we plot p s . Again, the B^ and Vt results are offset, because of the absence 
of photometry for bright stars. The abrupt decrease in the V band normalized dipole moment at 
the faint end of the Skymap catalog is probably incorrect. It is unlikely that the fainter stars and 
DGL would conspire to produce an anisotropy directed opposite to that resulting from the brighter 
stars; thus it appears that the anisotropy in B (V) is at least 5% (3%). For B, the (right ascension, 
declination) of the anisotropy remain quite constant for limiting magnitudes of 4 through 9, with 
value (115°, —43°). For V, the anisotropy direction gradually shifts to (115°, —51°) at magnitude 
10.5. 

The best available catalog for studying the anisotropy in the UV was constructed by the S2/68 
experiment (Boksenberg et al. 1973) on the European Space Research Organization's TD-1 satellite, 
which provided broadband flux measurements in four channels. The full catalog was not published, 
but Gondhalekar, Phillips, & Wilson (1980) presented tables of the integrated starlight in hundreds 
of patches, covering the entire sky. Gondhalekar et al. also estimated the flux contribution from 
the DGL and from stars too faint to be included in the catalog, using background data. In some 
cases, two or more stars were observed simultaneously. Such blended stars were excluded from 
the catalog, and Gondhalekar et al. estimated their contribution to the total flux, although the 
distribution of these unresolved stars on the sky was not reported. The Mathis et al. (1983) ISRF 
was constructed to maximize consistency with the Gondhalekar et al. (1980) results and with other, 
limited-coverage UV surveys. The adopted ISRF fluxes are somewhat higher than the total fluxes 
of Gondhalekar et al. 

In estimating the dipole in the UV, we use tables of integrated starlight in 10° x 10° patches in 
Gondhalekar (1989). These differ somewhat from the fluxes in Gondhalekar et al. (1980), because 
the absolute calibration of the S2/68 experiment was revised. In Table 3 we give vu v and the 
normalized dipole moment p s for the four UV bands; the listed dipole moments are likely lower 
limits. Note that the normalized dipole moment is substantially greater in the UV than in the 
visible. This difference presumably results from the differences in the stellar populations giving rise 
to the two spectral regions and the greater importance of extinction by dust in the UV. The UV 
radiation is dominated by early-type stars, which are relatively few in number and hence less evenly 
distributed than main sequence stars in general, all of which contribute to the visible radiation. In 
Table 4 we display the direction of the anisotropy for all of the considered bands. 
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For simplicity, we will assume that the anisotropy is independent of wavelength when calculat- 
ing drift speeds in the following section. We will adopt an anisotropy of 10%, intermediate between 
the anisotropics inferred for the visible and the UV. 

7. Gas-Grain Drift in the Diffuse ISM 

We will estimate drift speeds for grains in the cold neutral medium (CNM; T = 100 K, ran = 
30 cm -3 , n e = 0.045 cm -3 ) and the warm neutral medium (WNM; T = 6000 K, bh = 0.3 cm -3 , 
n e = 0.03 cm -3 ). 

In the presence of a magnetic field B, the steady state drift velocity is given by 
i \ F ( srn $ LOTs'mO \ 

(V\\,V±,V±±) = — cos6>, — — — r2,— 7 — ^ ( 37 ) 

where F is the net force excluding the magnetic and drag forces, Fd ra g = — ^Vd, 6 is the angle 
between B and F, lot = ZeB/cn is the product of the gyrofrequency and gas-drag time, and 
the subscripts (||, _L, _L _L) denote components (parallel to B, perpendicular to B and in the plane 
spanned by B and F, perpendicular to the plane spanned by B and F). Thus, drift across magnetic 
field lines is suppressed by a factor [1 + (lot) 2 ] 1 / 2 . 

We display F pc and F pc j for the CNM and WNM in Figure 15, and -Fdrag.H and ^drag,Coui 
in Figure 16. In order to estimate the largest effect that photodesorption is likely to have, we 
adopt 6> ra d = 90° (i.e. spin axis perpendicular to the net radiation flux); in this case F p d || S and 
Fdrag, H || — S. In order to gauge the importance of the photoelectric and photodesorption forces, we 
compute drift speeds both with and without including them. Specifically, we consider the following 
three cases: 

Case I .F ra d = F drag H + Fd rag ,He + Fd ragi coui , (38) 
Case II F rad + F pc = F° rag H + Fd rag) He + -^drag.Coui , (39) 

and 

Case III F rad + F pe + F pd = Fd rag) H + F dragi Hc + ^drag,Coui • (40) 

The resulting drift speeds are plotted in Figures 17 and 18. For a > 0.1/um carbonaceous grains, we 
find drift speeds v,} > 0.01 kms -1 in the CNM, and > 0.15kms -1 in the WNM; for a > O.l^m 
silicate grains, we find Vd > 0.002 kms -1 in the CNM, and Vd > 0.1 kms -1 in the WNM. The 
softness of the ISRF prevents F pe and F p d from substantially increasing the drift speeds over those 
that would hold if only F ra d were acting. These drift speeds were determined under the assumption 
that the radiation flux S is along the magnetic field. As seen in Figure 19, the suppresion factor 
for drift perpendicular to the field is substantial even for micron-sized grains. 
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7.1. Consequences 

Since grains of different sizes drift with different speeds, grains collide, with coagulation likely 
at low collision speeds. As a — > 0, v& — > 0, so that if grains always stick upon colliding, the timescale 
t a for a very small grain to become attached to another grain is given by 

t^WTr / a 2 u d (a)— ^da. (41) 

Here we estimate the maximum likely values for t a in the CNM and WNM, by adopting the drift 
speeds computed for case III above (i.e., including surface processes and assuming S || B; t a is 
only slightly longer for cases I or II). We adopt the grain size distribution from Weingartner & 
Draine (2001a) with Ry = 3.1 and be = 6 x 10~ 5 (their favored distribution for the average diffuse 
ISM). For coagulation with the carbonaceous grain population, t a ~ 8 x 10 9 yr (6 x 10 10 yr) in the 
CNM (WNM). For coagulation with the silicate population, t a ~ 1.5 x 10 10 yr (4 x 10 10 yr) in the 
CNM (WNM). Thus, drift due to forces associated with anisotropic radiation apparently results in 
negligible coagulation in the diffuse ISM. 

Gas-grain drift can also lead to variations in the dust-to-gas ratio. If the size of a cold neutral 
medium region is ~ lpc, then the time to separate micron-sized grains from the region could be as 
short as « 2 x 10 7 yr, which is perhaps comparable to the cloud lifetime. 



8. Summary 

We have investigated two forces that act on dust grains exposed to anisotropic radiation 
fields, due to the asymmetric photon-stimulated emission of particles from the grains. First, we 
have improved on Lafon's (1990) treatment of the photoelectric force, due to the recoil associated 
with photoelectrons, by using the realistic prescription of Kerker & Wang (1982) for estimating 
the emission asymmetry and employing the recent photoemission model of Weingartner & Draine 
(2001b). 

The photoelectric force depends on the ambient conditions relevant to the grain charging, 
including the shape of the radiation spectrum and the parameter GVT/n e . We find that the 
photoelectric force is comparable to the radiation pressure in regions where the grain potential 
is relatively low and the radiation spectrum is relatively hard (Figures 4-6). We find that the 
photoelectric force is unimportant compared with the radiation pressure in the cold and warm 
neutral media for grains larger than a few hundred A, because of the softness of the interstellar 
radiation field, as estimated by Mezger et al. (1982) and Mathis et al. (1983). 

We have also investigated the force resulting from the asymmetric photodesorption of H atoms 
adsorbed on the grain surface. The surface physics and chemistry of likely grain materials are not 
currently well-known. We have adopted a very simple, yet plausible, model in order to examine the 
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potential importance of this force. In particular, we assume that H chemisorbs at all surface sites 
with sufficient binding energy to supress thermal desorption and site-to-site diffusion. 

In our model, the total force associated with gas-phase atoms which collide with the surface 
includes contributions from three processes: photodesorption, H2 formation, and elastic reflection 
of atoms incident on the grain surface. We call the force on a grain that is not drifting with respect 
to the gas the "photodesorption" force. When the grain drifts, the total force (due to the above 
three processes) minus the photodesorption force equals the drag force, which is thus modified from 
the traditional treatment considering only the reflection of incident atoms (§5). 

With our standard choices for parameter values (Table 2), we find that the photodesorption 
force dominates the radiation pressure for grains with a > 0.1/zm when the radiation field is 
relatively hard (Figure 9). Laboratory studies of the surfaces of candidate grain materials, including 
determinations of photodesorption cross sections, are needed in order to clarify the true importance 
of this force. 

Grains exposed to anisotropic radiation will drift relative to the gas due to the combined 
action of the radiation pressure, photoelectric force, and photodesorption force. We find that the 
radiation field in the diffuse ISM is ~ 10% anisotropic in the visible and ultraviolet (§6). In §7, 
we estimate the maximum likely drift speeds for grains in the cold and warm neutral media. We 
find the timescale for grain coagulation resulting from this drift to be > 5 x 10 9 yrs. Large grains 
(a ~ 1/xm) might drift rapidly enough to be substantially separated from the gas in cold diffuse 
clouds. In a separate paper (Weingartner & Draine 2001c), we will investigate gas-grain drift in 
photodissociation regions. 

This research was supported in part by NSF grant AST-9619429 and by NSF Graduate Re- 
search and International Fellowships to JCW. We are grateful to R. H. Lupton for the availability 
of the SM plotting package. 
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Fig. 1. — Radiation pressure efficiency factors for neutral carbonaceous and silicate grains, averaged 
over the interstellar radiation field (ISRF) and blackbody spectra with indicated color temperatures 
(cut off at 13.6 eV). The kink at a = 50 A results from the Li & Draine (2001) prescription for 
blending PAH and graphite optical properties. 
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Fig. 2. — The asymmetry factor A as a function of the incident photon energy hv, for carbonaceous 
grains. The grain radius is indicated for each curve. 
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Fig. 3. — Same as Figure 2, but for silicate grains. 
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Fig. 4. — Fpe/Fj-ad for carbonaceous grains, T c = 3 x 10 4 K, and two gas temperatures: 100 K (solid) 
and 1000 K (dashed). The values of G^/T/n ej in cm 3 , are indicated. Curves labelled CNM 
(WNM) are for cold (warm) neutral media and the ISRF. 
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Fig. 5. — Same as Figure 4, but for silicate grains. 
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Fig. 6. — F pe /F ra d for carbonaceous grains, T = 100 K, GVT/n e = lO 3 ^/ 2 cm , and four values 
of T c as indicated. In all cases, the radiation field is assumed to be cut off at 13.6 eV . 
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Fig. 7. — Normalized rate for photodesorption for carbonaceous grains, a blackbody radiation 
spectrum (cut off at 13.6 eV) with T c = 3 x 10 4 K, and several values of grain size a, as labeled. 
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Same as fig. 7, but for silicate. 
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Fig. 9. — Lower panel: The photodesorption asymmetry integral appearing in the expression for 
the maximum possible force for flipping grains, (eq. 23). For carbonaceous (silicate) grains, 
the curve corresponding to a blackbody radiation field with T c = 3 x 10 4 K (cut off at 13.6 eV ) lies 
below (above) that for the ISRF. Upper panel: Ratio of the maximum possible photodesorption 
force to the radiation pressure force, for parameter choices as indicated in Table 2 (note that Eu 2 
is not relevant here). 
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Fig. 10. — Upper panel: The photodesorption force as a function of angle # ra d between Co and S, 
normalized to [see eq. (23)]. Lower panel: 9{ is the angle between F p d and S. 
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Fig. 11. — The occupation fraction / as a function of the angle 6 with respect to the spin axis, for 
carbonaceous grains with a = 0.3/um, sh = s c = 1, pu 2 = Pu, a dilute blackbody radiation field 
with T c = 3 x 10 4 K, and two values of a: a = 0.01 and a = 1. 
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Fig. 12. — Fpd/F^ and cos#f for non-flipping carbonaceous grains with a = 0.3^m, sh = 1, 
s c = 0.1, pu 2 = PHj a blackbody radiation field with T c = 3 x 10 4 K, and various values of (a, 
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mag 

Fig. 13. — Energy density per wavelength interval, normalized to that of the ISRF, as a function 
of limiting magnitude. Results are from the Skymap and Tycho star catalogs. The Tycho fluxes 
are systematically low due to the absence of photometry for thousands of bright stars. 
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Fig. 14. — The normalized dipole moment p s , i.e., the fraction of the energy density in the ISRF 
that is in the anisotropic component, as a function of limiting magnitude. The Tycho anisotropies 
are systematically low due to the absence of photometry for many of the brightest stars. 



Table 1. Momentum Transfer Rates, Per Site 



process rate 

photo desorption |i? p dPH/ 
H2 formation |-Rarr s cPH 2 / 



reflections [1 — su + /(sh — s c )] J °° dvP(v)nul 2 v(2m p v) = [1 — sh + f{su — s c )] m p R\ 
sticking [s H - f(su - s c )] J °° dvP(v)nYtl 2 v(m p v) = [s H - /(sh - s c )] mp f arr 







Table 2. Canonical Photodesorption Parameters 



parameter 


symbol 


value 


photodesorption rate (graphite) 


R l& 


5 x lO-^Gs" 1 


photodesorption rate (silicate) 


R° 
n P d 


2 x lO-^Gs -1 


kinetic energy of photodesorbed H 


En 


2eV 


kinetic energy of newly-formed H2 




leV 


surface site area 


f 


6 A 2 


sticking probability 




1 


combination probability 


Sc 


1 
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Fig. 15. — The photoelectric and photodesorption forces for grains in the warm and cold neutral 
media. 
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Fig. 16. — The drag forces due to collisions with neutral H and to Coulomb interactions with 
distant ions. 




Fig. 17.— 



The drift speed for carbonaceous grains in the warm and cold neutral media. 
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Fig. 18. — The drift speed for silicate grains in the warm and cold neutral media. Coincidentally, 
the minimum value of -F ra d + F pc + F p d in the CNM is very close to zero; hence the sharp drop in 
Vd when a « 350 A. 
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Fig. 19. — The product lot, where to is the gyrofrequency for the charged grains in the ambient 
magnetic field, and r is the drag time. Drift transverse to the magnetic field is suppressed by a 
factor [1 + (lot) 2 ] 1 / 2 [see eq. (37)]. 
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Table 3. UV Energy Densities and Dipole Moments 



Quantity 


2740 A 


2365 A 


1965 A 


1565 A 


starlight dipole {p s vu y ) 


& ' b 1.27 


0.843 


1.30 


1.67 


total starlight (vu v )* 


a,6 443 


3.23 


4.73 


5.27 


total vu v a,c 


4.90 


5.47 


6.47 


6.90 


ISRF vu v a ' d 


8.30 


6.20 


7.00 


8.20 


dipole/ISRF 


0.152 


0.136 


0.185 


0.205 



a 10 14 erg cm 3 
b Gondhalekar 1989 
c Gondhalekar et al. 1980 
d Mathis et al. 1983 



Table 4. Anisotropy Directions 



Band 


RA 


dec 


/ 


b 


1565 A 


141° 


-62° 


281° 


-8.2° 


1965 A 


132° 


-58° 


275° 


-8.9° 


2365 A 


124° 


-55° 


270° 


-10.5° 


2740 A 


129° 


-56° 


273° 


-9.3° 


B 


115° 


-43° 


257° 


-9.8° 


V 


115° 


-51° 


264° 


-13.6° 



Note. — Celestial coordinates: right ascension (RA) and declination (dec); Galactic coordinates: 
longitude (/) and latitude (b). 



